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In memoriam of Martine Babillot. 

Abstract. We prove Gaussian tail estimates for the transition probability of 
n particles evolving as symmetric exclusion processes on "L d , improving results 
obtained in We derive from this result a non-equilibrium Boltzmann-Gibbs 
principle for the symmetric simple exclusion process in dimension 1 starting 
from a product measure with slowly varying parameter. 



1. Introduction 

To derive sharp bounds on the rate of convergence to equilibrium is one of 
the main questions in the theory of Markov processes. In the last decade, this 
problem has attracted many attention in the context of conservative interacting 
particle systems in infinite volume. Fine estimates of the spectral gap of reversible 
generators restricted to finite cubes and logarithmic Sobolev inequalities have been 
obtained. We refer to [H] for the recent literature on the subject. From these bounds 
on the ergodic constants, polynomial decay to equilibrium in L 2 has been proved 
for some processes. For instance, Bertini and Zegarlinski £Q, |2] proved that the 
symmetric simple exclusion process in 1 d converges to equilibrium in L 2 at rate 
t~ d l 2 . For a class of functions / that includes the cylinder functions, there exists 
V(f) finite such that 

\\Ptf- < / >« || 2 < 1 (i i 



12 - {i + tyi 2 

for all t > 0. Here Pt stands for the semi-group, < / > a for the expectation of 
/ with respect to the Bernoulli product measure with density a and ||/||2 for the 
L 2 norm of /. Janvresse, Landim, Quastel and Yau [HJ and Landim and Yau |10| 
extended the algebraic decay in L 2 for zero range and Ginzburg-Landau dynamics. 

We refine in this article the Gaussian upper bounds obtained in [Hj for the tran- 
sition probabilities of finite symmetric simple exclusion processes evolving on the 
lattice Z d . Our approach is based on a logarithmic Sobolev inequality and on Davies 
0] method to derive estimates for heat kernels. 

Consider n > 2 indistinguishable particles moving on the d-dimcnsional lattice 
Jj d as symmetric random walks with an exclusion rule which prevents more than 
one particle per site. The dynamics can be informally described as follows. There 
are initially n particles on n distinct sites of Z d . Each particle waits a mean one 
exponential time at the end of which, being at x, it chooses a site y with probability 
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p(y — x), for some finite range, irreducible, symmetric transition probability p{-)- 
If the site is vacant, the particle jumps, otherwise it stays where it is and waits a 
new mean one exponential time. 

The state space of this Markov process, denoted by £ n , is the collection of all 
subsets A of Z d with cardinality n: 

£ n = {AcZ d , \A\ = n} ; 

while its generator C n is given by 

(C n f)(A) = £ p(y-x)[f(A XtV )-f(A)}, (1.1) 

where A x y stands for the set A with sites x, y exchanged: 

f (A\{x})U{y} ifxeA,y?A, 
A x , y = I (A \ {y}) U{x} if y e A, x <? A, 
[ A otherwise. 

In formula (|l.l|l summation is carried over all bonds {y, z} to avoid counting twice 
the contribution of the same jump. 

It is easy to check that the counting measure on £ n , denoted by v n (v n (A) = 1 
for every A in £„), is an invariant, reversible measure for the process. 

Fix a set A in £ n and denote by f t = ff° the solution of the forward equation 
with initial data Sa ■ 

( d t ft = Cnft , 

\ f (A) = 1{A = A } . { ■ > 

The main result of the article provides a Gaussian estimate for the transition 
probability f t . Denote by x = (xj., . . . , x n ) the sites of (Z d )™. For a configuration 
x, let Xjj be the j-th coordinate of the i-th point of x: Xj ; j = Xi -ej, where • stands 
for the inner product in M. d and {ei, . . . , e<j} for the canonical basis of M. d . The 

2 ^2 



Euclidean norm of (R ) n is denoted by ||x|| so that ||x|| = 



Denote by $ the Legendre transform of the convex function w 2 coshiu: 
&(u) = sup {uw — w 2 coshw} . 

An elementary computations shows that Q(w) ~ w 2 for w small and $(w) ~ w \ogw 
for w large. 

Theorem 1.1. Fix a set Aq — {zi, . . . , z n } in £ n . Let f t be the solution of the 
forward equation jl.ty) . There exist finite constants C'2 — C2(n,d,p), ao = ao(p) 
such that 

Jt{A) _ 2_, (1 + T) nd/2 CX P1 2 (logT) 2 V a 2 T 

for every T > C2 and every set A = {xi, . . . ,x n }- In this formula, summations 
is performed over all permutations a of 11 points and x CT stands for the vector 

( x cr(l) : ■ ■ ■ : x a{n))- 

The asymptotic behavior of $(■) at the origin shows that for every 7 > 0, there 
exists a constant a\ = ai(p, 7) such that 
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for all T > C2 and all sets A such that ||x CT — z|| < 'yT/logT for all permutations 
a. Furthermore, since 

/ ||x CT -z||logT \ lA /||x CT(i) -Zi||logT 

V nlT I — „.2-*i 



aiT ) 11 V o3,T 
i=i 

we have that 

C 2 _/ a T A^^llx^) - Zi||logT> 

i=l 



, . >r^ ^2 r a T ^ /||x . w -z i ||logTN-i 



For a fixed 7 > 0, in last formula we may replace by C(j)w 2 if Hx^) — Zj|| < 

a^Tj/logT and by C(7)w log it; otherwise. 



2. BOLTZMANN-GlBBS PRINCIPLE 

We prove in this section the Boltzmann-Gibbs principle for the symmetric sim- 
ple exclusion process out of equilibrium in dimension 1. This result allows the 
replacement of average of local functions by functions of the empirical density in 
the fluctuation regime and is the main point in the proof of a central limit theorem 
around the hydrodynamical limit for interacting particle systems (cf. 6 , Chap. 
11). We restricted our attention to dimension 1 because Lemma IT^I below has only 
been proved in d = 1 . 

The Boltzmann-Gibbs principle for one-dimensional processes out of equilibrium 
was proved in [2] through the logarithmic Sobolev inequality. A different version, 
involving microscopic time integrals, is presented in [H| and uses sharp estimates on 
the comparison between independent random walks and the symmetric exclusion 
process. 

Fix a profile po :R d ^ [0, 1] in C 1 (R d ) with a bounded derivative and consider a 
sequence of product measures {y^i.y N > 1} on {0, 1} Z<1 associated to this profile 
so that 

E vN Mx)\ = Po (x/N) . 

po() 

Denote by P„n the probability measure on the path space D(K + , {0, 1} Z ) corre- 

PO < ■ ) 

sponding to the symmetric simple exclusion process starting from ^(.) and speeded 
up by N 2 . Expectation with respect to P„n is denoted by . 

po() PO(') 

For x in Z d , let 

p N (t,x) = E„n Mx)] . 

po() 

Of course, p N (t,x) is the solution of the linear equation 

' d tP N (t,x) = n 2 Y, p(v - [A*, y) - /Am)] , 

< y ez d 
^p N (0,x) = p Q (x/N) . 

This equation can be written as d t p N (t, x) — N 2 dp N (t, x), where L\ is the gener- 
ator introduced in l|l.l|) ■ Next proposition is the main result of this section. 
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Proposition 2.1. Let d = 1 and fix T > 0, a finite subset AofL such that \ A\ > 2 
and a continuous function H in L 1 (M). Then, 



lira E„w 

N^oc po(-) 



( [ * ^tt2 e # o*/^ n + z ) - ^(*> x + z )i 



2n 

= 



The proof of this proposition is presented at the end of this section. The 
Boltzmann-Gibbs principle is a simple consequence but requires some extra no- 
tation. 

For a finite subset A of Z and < a < 1 , let 



${A,a) = Y[[ri(x)-a] 



xeA 



By convention, we set ^(cf>,a) = 1. Each cylinder function / : {0, 1} Z — > R can be 
written as 



f(r,) = £f(A,a)*(Aa) . 

Ae£ 

A straightforward computation shows that for each finite set A, f(A, •) is a smooth 
function, in fact a polynomial. 

For a cylinder function /, let / : [0, 1] — ► K be the real function defined by 
/(a) = E Va [f(r])} and let 

I/fa, a) - /(ry) - /» - /'(a)fo(O) - a] . 

Note that f(</>,a) = /(a) and that Szez Ki^L a ) = f'( a )- I n particular, it follows 
from a simple computation that 

r f ( v ,a) = £ f(M» -*»(°)] + E E f(Aa)*(A«)- 

zGZ n>2A6£„ 

Fix a smooth functions iJ in L 1 (R). By the previous formula, 
*J2H(x/N)r f (T x r),p N (t,x)) 

z<EZ 

= E ^7172 E (*/tf)f({*}. p w (^ + *) - v(x)\ 

z6Z xGZ 

+ E E j^^H(x/N)ftA,fF(t,x))T w ilf(A,fr(t,x)). 

n>2A<E£ n zSZ 

Note that the sums in 2, n and A are finite because / is a cylinder function. Since 
H, p N (t, •) and f({z}, •) are smooth functions, a change of variables shows that the 
first term is of order TV -1 / 2 and that the second is equal to 

E E ^Ew)p,/(M))n^+ z )-A^+^)] + o(^v- i/2 ), 

n>2Ae£„ zeZ z£A 

which is exactly the expression appearing in Proposition 12.11 Since f(A, •) and 
p N (t, •) are smooth bounded functions, we have proved the following result, known 
as the Boltzmann-Gibbs principle. 
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lim E„w 

N~>oo po(-) 



Corollary 2.2. Let d = 1 and fix T > 0, a cylinder function f and a smooth 
function H in L (K). Then, 

^ dt^Y, H ^l N ) T f^P N{t ^))) 2 ] = °" 

The proof of Proposition 12. II is based on three lemmas concerning the decay of 
the space-time correlations of the symmetric exclusion process. We start with a 
general result which will be used repeatedly. 

Fix n > 2 and denote by f t (A, B) = f t N (A, B) the semi-group associated to the 
generator N 2 C n . For a finite subset A of Z, let 

1(A) = ^2p(v-x). 

x,y£A 

Note that 1(A) vanishes unless A contains two sites which are within a distance 
smaller than the range of the transition probability. Next lemma follows from 
Theorem 11.11 and a straightforward computation. 

Lemma 2.3. For all T < oo, n > 2, there exists a finite constant C3, depending 
only on n, p and T such that 

jf^ d + J^m/a £ f t N s (AB)I(B) < C 3 A N (m,t) 

for all A in E n , N > 1, < t < T ; where 

{N- 1 ifm = 0, 

N~ 2 if m=\, 

log N/N 2 Vl + tN^ ifm = 2, 

l/N 2 VT+tN 2 ~ ifm>3. 

We now introduce the space-time correlations, also called w-functions in [7]- For 
a finite subset A of Z and t > 0, let <p N (t, <f>) = 1, 

v N (t,A) = E v n [ JJ {ijt( x ) - P N (t> x )} • 
xeA 

Notice that ip N (t, {x}) vanishes for all x. An elementary computation shows that 
f dw N (t,A) = N 2 C nV N (t,A) + G N (t,A) , 
\^(0,A) = 0, (2 '" 
where n = \A\ and G N (t, A) is given by 
N 2 J2 p(z-y){<p N (t,A\{z})-^ N (t,A\{y})}{ P N (t,{z})-p N (t,{y})} 

y,zeA 

-(N 2 /2) J2 P (z-y) V N (t,A\{y,z}){p N (t,{z})- p N (t,{y})) 2 . 

y.zeA 

Here again summation is carried over all bonds. Notice that the first line vanishes 
for n = 2 and that the second line vanishes for n = 3. 

The linear differential equation 1|2.1|) has a unique solution which can be repre- 
sented as 

ip N (t,A) 



[ ds V / t _ s (4,B)G"( S ,B) 
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so that the space-time correlations ip N (t,A) can be estimated inductively in n. 

Next lemma is due to Ferrari, Presutti, Scacciatelli and Vares 0. In the proof 
of Proposition 12 . 1 1 we don't need such sharp estimates. 

Lemma 2.4. Assume that d = 1 and fix T > 0. For each n > 1, there exists a 
finite constant C4 = Ci(n,p, po,T) such that 

sup \<p»(t,A)\ sup \^(t,A)\ < ^g^. 

0<t<T 0<t<T 

For < s < t < T, 1 < k < n and B G £ k , A e £ n , let 

R N (s,A;t,B) = E v » () [l[[T, a (x)-p N ( S ,x)] ]\h(y) - P N (t,y)] 

xeA yeB 

Since s and A will be fixed, most of the time, we denote Rn(s, A; t, B) by R,N(t, B). 
Notice that in this definition we do not require A and B to have the same cardinal- 
ity. An elementary computation shows that Rjy(t,B) is the solution of the linear 
differential equation 

( d t R N (t, B) = N 2 C k R N {t, B) + H N (t, B) for t > s , 
\R N (s,B) = J N (s,A,B) , 

where k = \B\, 

J N {s,A,B) = E <o{) [l[[n s (x) - p N ( S ,x)}Y[[ Vs (y) - p N (s,y)} 

xeA yGB 

and Hff(t, B) is given by 
N 2 p(z-y){RN(t,B\{z})-R N (t,B\{y})}{p N (t,{z})-p N (t,{y})} 

y,zeB 

-{N 2 /2) Y p(z~y)RN(t,B\{y,z}){p N (t,{z})-p N (t,{y})} 2 . 

y,zeB 

Notice that i?jv(t, </>) — p N (s,A), that H^(t,B) vanishes for n — 1 and that 
Jn(s, A, B) is not equal to tp N (s, AU B) but given by 

Yf N (s,(AAB)UC)l[[l-2p N (s,{x})] J] p N (8,{x})[l-p tf (s,{x})], 
c xec xe(AnB)\c 

where the summation is carried over all subsets C of A n B and where AAB stands 
for the symmetric difference of A and B. 

The differential equation 12.21) has a unique solution which can be represented 

as 

R N (t,B) = Y fr(B,C)J N (s,A,C) + J du Y fr-u(B,C)H N (s + u,C) , 
where r = t — s. This last notation is systematically used below. Let 



U N (t,B) = Y fr(B,C)J N (s,A,C) 



ce£ k 
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Lemma 2.5. Fix 2 < k < n, < s < t < T and A in S n ■ There exists a finite 
constant C4 = C 4 (p, n, T, po) such that 

C 4 B(n - k) 



sup \U N (t,B)\ < 



Be£ k 



(l + rN 2 ) k / 2 



where B(2j) = and B(2j + 1) = log N/N^> +1 for j > 0. 
Proof. UN(t,B) is absolutely bounded by 



J2j(s,j) fr{B,D U E) , (2.4) 



j=0 DcA.DeSj 



where 



J(s,£) = sup \Jn(s,A 7 C)\ < max sup lip (s,D)\ , 

C££ k m D£E m 

\Cr\A\=l 

where the maximum is carried over n+k— 21 < m < n+k+£. Last inequality follows 
from the explicit formula for Jpj(s,A, C). By Lemma |2 .41 the previous expression 
is less than or equal to C 4 B(n + k — 2£). On the other hand, by Theorem ll.il 



D<zA,D££j 



Therefore, 

1/7 (t ml < r ^ B(n + k-2j) C 4 B(n-k) 
\ U ^ B )\ $ ( 1 + rN 2y/2 S {1 + rN 2 )k /2 ■ 

This concludes the proof of the lemma. □ 

We are now in a position to prove the main result towards the Boltzmann-Gibbs 
principle. 

Lemma 2.6. Fix n > 2, there exists a finite constant C4 = C 4 (n,p,T, po) such 
that 

logiV 1 -1 



\R N (s,A;t,B)\ < C 4 { ; 



N 2 1 + (t - s)N 2 J 
for all A, B in £ n , < s <t <T. 

Proof. Fix n > 2, s > and A in £ n . For 1 < k < n, denote by RN,k(t,-) 
the solution of the linear differential equation (|2.2|l . Since the equation for Rjsr,k 
involves R^.k-i, RN,k-2, an induction argument on k is required. A simple pattern 
appears only for k > 7. Hence, for 1 < k < 6, we need to proceed by inspection, 
making the proof long and tedious. 

Consider k = 1. In this case Hn vanishes and, by Lemma 12.51 

\r n &{x})\ = \u N ( t ,{x})\ < {1+ ; N2)1 ; 2 ■ 

Here and below {dj,j > 1} are finite constants depending on n, p, T and po which 
may change from line to line. 
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For k = 2, since i?jv(t, 4>) is time independent and absolutely bounded by B(n), 
the previous estimates and Lemma 12.51 show that 

NB(n — 1) 

rN 2 ) 1 / 2 

Therefore, by the explicit formula l|2.3[) for R^^itjB) and by Lemmas 12.31 12.51 

(B{n-2) B(n-l) 



f NB(n — 1) l 
\B»&B)\ < ^{ (1 + l N J /2 + B(n)}l(B) 



Ir u ml <- f Wji-X) B n-i -l 



because S(n) < B(n — 1). Notice that this inequality proves the lemma for n = 2 
because Rn,2^^ B) = Rn(s, A; t, B). 

The estimates for Rn,i and Rn,2 give bounds for Hn. 3 which in turn, together 
with the explicit formula 1)2. 3|) for R]y^(t,B) and Lemmas 12 . 31 12 . 51 show that 

In r* ml <- f gtn^j) i?(»-2)logjV 1 
\R N , 3 (t,B)\ < a 3 ( (1+riv2) 3 /2 + A/(1 + rAr2)1/2 | 

Here we used the fact that £?(n — 3) = NB(n — 1) to eliminate one of the terms 
appearing in the expression of RN^(t, B). 

We repeat this procedure for k — 4, 5 and 6 to obtain that 

I I r B{n-A) B{n-3) B{n-2)logN\ 
\R NA (t,B)\ < a^ {i + rN2)2 + N{l + rN2)1/2 + W2 ) , 

IB U ml < f B(n-5) B(n - 4) B(n-3) \ 

\Hn, 5 ^B)\ < «5( (1 + rAr2)5/2 + Ar(1 + rJv2)1 / 2 + N 2 ) > 

IB ^ ml <r f g ("~6) g(n-5) g(n-4) ^ 

1^6^5)1 < «6( (1 + r7v2)3 + N{l + rN2 y/2 + -^^) ■ 

For k = 5, we used the fact that log N < B(£ - 1). 

A pattern has been found for k = 5, 6. It is now a simple matter to prove by 
induction that this pattern is conserved so that 

B(n-k) B(n-k + l) , B(n~k + 2)\ 



\RN,k{t,B)\ < a 6 {— - 



(l + rN 2 ) k / 2 Nil + rN 2 ) 1 / 2 N 2 J 

for k > 7. It remains to recall the definition of B(j) and to recollect all previous 
estimates to conclude the proof of the lemma. □ 

Notice that we could have set B(l) = N^ 1 for the estimates in the previous 
lemma. Taking B(l) = log N/N simplifies slightly the notation since we have that 
B(n + 2) = B^N^ 1 for all n > and we miss only a logiV factor, which is 
irrelevant for our purposes. 

We are now in a position to prove Proposition ^. II With the notation introduced 
in this section, the expectation appearing in the statement of the proposition be- 
comes 

1 H(x/N)H(y/N) f dt f ds R N (s, A + x;t, A + y) , 

where A + x is the set {z + x : z € A}. By Lemma \2. 61 and a change of variables, 
this expression is bounded above by 

N \N ^ 1 v 1 ' 
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which proves Proposition l2.il 

We conclude this section with an observation. The same arguments presented 
above in the proof of Proposition 12 . II shows that 



E, 



poCO 



< 



+ 



2 f dt f ds^-Y, F{ ~P N ^ x )) H ^l N )^sH){x/N) 

J ° J ° iV x 6 Z 

C{ PlPOl T)(^Y.\ H{ - x l N )\ 



for some finite constant C(p, po, T). Here F(a) = a(l — a). 



3. Gaussian tail estimates for labeled exclusion processes 

Fix n > 2 and a finite range, symmetric and irreducible transition probability p(-) 
on Z d . Consider n labeled particles moving on the <i-dimensional lattice Z d through 
stirring. This dynamics can be informally described as follows. The n particles start 
from n distinct sites of Z d . For each pair (x, y) of Z d , at rate p(y — x), particles at 
x, y exchange their positions. This means that if there is a particle at x (resp. y) 
and no particle at y (resp. x), the particle jumps from x to y (resp. from y to x). 
If both sites are occupied, the particles change their position and if none of them 
are occupied, nothing happens. 

The state space of this Markov process, denoted by B n , consists of all vectors 
x = (sci, . . . ,x n ) of (1 d ) n with distinct coordinates: 

B„ = jx = (xi, . . . , x n ) £ (Z d ) n : Xi ^ xj for i ^ jX 

while the generator L n is given by 

(i»/)(x) = ^ p(y-xMo*"x)-f(x)]. (3.1) 

In this formula, for a configuration x = (x\, . . . , x n ) in B n , a x ' y x is the configura- 
tion defined by 

{y if xi = x, 
x if Xi = y, 
Xi otherwise. 

This generator corresponds to the generator p. 1(1 in which particles have been 
labeled and are therefore distinguishable. 

It is easy to check that the counting measure on B n , denoted by /x n , is an 
invariant reversible measure for the process. The goal of this section is to obtain 
sharp estimates on the transition probability of this Markov process. To state the 
main results of the section, fix a state z in B n and denote by f t the solution of the 
forward equation: 

J d t ft = L* n f t , , . 

\ / (x) = l{x = z}. W 

Recall that we denote by $ the Legendre transform of the convex function 
w 2 coshw. 
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Theorem 3.1. Fix n > 1 and a point z = (z\, . . . , z n ) in B n . Let ft be a solution of 
the forward equation 13.2(1 . There exist finite constants Ci = C 2 (n,d,p), ao — ao(j>) 
such that 

/or every T > C 2 and every configuration x. 

Since $(u>) ~ uj 2 for w small, for 7 > 0, there exists a finite constant a\ = 
a\(jp, 7) such that 

M X ) < ( l +T )nd/2 CX P{ fll r / 

for every T > C 2 and every configuration x such that ||x — z|| < 7T/ logT. 
On the other hand, since x 2 cosha; < 2e 2x , <E>(u) > (u/2) log(it/4e). Hence, 

. , x ^ g 2 r ||x-z|| l|x-z||io g r 1 

/Hx) < {1 + T)nd/2 exp{ - ^-j^ log 4eflgr } 

for every T > C 2 . Of course this estimate is only interesting if |jx — z|j ^S> T / logT, 
Since the proof of Theorem 13.11 follows closely the one of Theorem 2.2 in 

we present only the main differences. Throughout this section, Co stands for a 

universal constant, which may change from line to line. 

We first need a logarithmic Sobolev inequality for the process Xt restricted to 

cubes. Fix an integer I and decompose the lattice 1 d into disjoint cubes {Afc : k > 

1} of length t 

Afc = Xk + {1, . . . , (} d for some x fe in Z d ; 
A fe n Aj ; = 4> for k^j and (J A fc = Z d . 

k>l 

For a vector k = (fc l7 . . . , fe n ), let Ak be the finite cube of (Z d ) n defined by Ak = 
Afc! x • • • x Afe n and let L\ k be the generator L n introduced in (|3.1|) restricted to 
the cube Ak. This means that jumps from Ak to its complement are forbidden as 
well as jumps from the complement to Ak- 

Lemma 3.2. There exists a finite constant C\ depending only on the transition 
probability p(-), the dimension d and the total number of particles n such that 

/(x)io g /(x) < c,e Yl {JW)-JW)} 2 1 (3-3) 

x£A k x,y£A k 

for all densities f with respect to the uniform probability measure over Ak- In this 
formula, the sum on the right hand side of the inequality is carried over all pairs 
x, y in Ak such that y = <t i,!, x for some x, y with p(y — x) > 0. 

Proof. It is well known that a symmetric random walk evolving on a d-dimensional 
cube satisfies a logarithmic Sobolev inequality of type l|3.3|l and that the super- 
position of independent processes satisfying logarithmic Sobolev inequalities also 
satisfies a logarithmic Sobolev inequality, the constant being the maximum of the 
individual constants. This proves l|3.3|) in the case where the cubes Afc are all 
different: ki ^ kj for i =/= j. 

It remains to consider the case where some cubes are equal. In this situation, the 
diagonal is forbidden because two particles cannot occupy the same site, and two 
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particles may exchange their position. Fix 2 < m < n and consider the hypercubc 
Ak = Afe x • • • x Afe of (Z d ) m . If we do not distinguish particles, we retrieve the 
symmetric simple exclusion process on A^ with m particles. This process satisfies a 
logarithmic Sobolev inequality of type l|3.3|) with a constant Co depending only on 
the dimension d and the transition probability p(-) jllj . It is not difficult to recover 
(I3.3|l for the random walk X t on Ak from this estimate. 

Indeed, let 2A fc , m be the subsets of A& with m points: £A fc ,m = {A C A& : 
\A\ = m}, let fJ-A k ,m be the uniform probability measure on EA fc ,m and, for a density 
/ : Ak — > M+ with respect to the uniform measure over Ak, let / : £A fc , m — > R+ be 
the density with respect to [iA k .m defined by 

/({iCi, . . . ,X m }) = —:'Y)f(x a ti),...,X a t m \), 

ml * — ' 

a 

where the summation is performed over all permutations a of m elements. 
With this notation, we may rewrite the left hand side of (|3.3|) as 

E f( A )EjfT, l ° S M +m] E f(A)logf(A), (3.4) 
Aes Afc ,m xgaAA> J( A ) Aes Ak , m 

where the summation over x is carried over all points x = (x\, . . . , x m ) such that 
{xi,.. .,x m } = A. 

It is not difficult to prove a logarithmic Sobolev inequality for the permutation 
of m points. Let § m be the set of all permutations a of m points. Consider the 
Dirichlet form D§ m defined by 

D Sm (g) = E ^) - g(t)} 2 ■ 

tr,5-eS m 

There exists a finite constant Co such that 

for all densities g with respect to the uniform probability measure on § m . 

Since f(x)/f(A) is a density with respect to the uniform probability measure 
over the set of all permutations, the first term is bounded above by 

Co E E {v/TW-VTw} 2 (3-5) 

for some finite universal constant. It remains to connect each point x in A to each 
point y in A by a path x = zo, . . . , z r = y such that Zj+i = <j x,v Zj for some x, y 
with p(y — x) > to estimate the previous term by the right hand side of l|3.3|l . 
This can be done as follows. 

Assume first that d = 1. To explain the strategy in a simple way, we allow two 
particles to occupy the same site in the construction of the path. The modifications 
needed to respect the exclusion rule are straightforward. Fix x and y in a same 
set A. Since both points belong to the same set, there exists a permutation a of 
m points such that yi = x a ^ for 1 < i < m. The path {zj} connecting x to y is 
defined as follows. We start changing the first coordinate x\ of x, keeping all the 
other constants, moving from x = (xi, . . . , x m ) to wi = [y\ = a; CT (i), X2, ■ ■ ■ , x m ). 
Note that the last configuration has two particles occupying the same site. At 
this point, we change the coordinate as^m, moving from a new configuration W2, 
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which is obtained from x, by replacing x\ by x a n\ and av(i) by avam, where 
(j = a o a. We repeat this procedure. If the orbit of 1 for the permutation a 
is the all set {1, . . . , m}, this algorithm produces a path from x to y. Otherwise, 
after completing the orbit of 1 by the map a, we choose the smallest coordinate 
not belonging to the orbit of 1 and repeat the procedure. 
Denote by r xy the path just constructed. Notice that 

1. its length is bounded by mi and 

2. all coordinates but one of each site z in r xy belong to the set {xi, . . . , x m }. 
Therefore, by Schwarz inequality, (|3.5|) is bounded above by 

c o E E l r ^l E {v/iW-v/TIbiT} 2 

< ctm^JJ^-JWi)} 2 E E • 

b Aes A , m x.yeA 

The last sum in the first line is performed over all pairs 6 = (61,62) of consecutive 
sites in the path r xy , while the first sum in the second line is performed over all 
pairs 6 = (61,62) such that 62 = a x,x+y for some x, y in Z d such that p(y) > 0. 
Since all but one coordinate of each site in r xy belong to {£1, . . . ,x m }, for each 
fixed bond 6 = (61,62) there is at most m£ possible sets A which might use this 
bond. For each set A, there is at most ml end points and ml starting points for the 
path. The last sum is thus bounded by 

< c «(m!) 2 ]T{v7(b27- V7(b7)} 2 - 

b 

This concludes the proof of the estimate of the first term in (|3.4() in dimension 1. 

The proof in higher dimension is similar. The idea is to consider a configuration 
x as a point in Z md and repeat the previous algorithm, moving the first coordinate 
of the first particle, then moving the first coordinate of the c(l)-particle, until all 
first coordinates of all particles are modified. At this point, we change the second 
coordinate of the first particle and repeat the procedure. This method gives a path 
of length at most Coimd and whose sites have all but one of the md coordinates 
equal to the coordinates of x. These two properties permit to derive the estimates 
obtained in dimension 1, replacing m by md. This proves that the first term in 
(|3.4I) is bounded above by the the right hand side of (|3.3|) . 

We focus now on the second term of l|3.4|) . By the logarithmic Sobolev inequality 
for m exclusion particles in a cube A^, this expression is less than or equal to 

Cm! £ ]T p(y){y/M x , x+y )- \[J{A)} 2 

for some finite constant C depending only on p(-) and d. By Schwarz inequality, 
this expression is bounded by the right hand side of (|3.3|) . □ 

The second main ingredient in the proof of Theorem 13.11 is an estimate of the 
action of the generator L n on certain exponential functions. 

For a vector 9 = {6\, . . . , 6 n ), 0, in R d , denote by ip$ the function ipg : B n — > R 
defined by Ve( x ) = exp{6* • x}. Here, x • y represents the inner product in (Z d ) n . 
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An elementary computation shows that there exists a finite constant ao, depending 
only on the transition probability £>(•), such that 

-^(L^Xx) < R(9), £^(,-z){M^-l} 2 < R(9) (3.6) 
for all x in B n , where 

d n 

R(9) = a ( cos Mao^,i} - l) • (3.7) 

j=i i=i 

Next result relies mainly on Lemma 13.21 and on the bounds (|3.6(l . Its proof 
follows closely the one of Lemma 4.3 in [§] and is therefore omitted. For a positive 
function ip : B n — * R, denote by the Dirichlet formula defined by 

S)^(u) = (1/2) J] ^ p(y - ^{ttK^x) - ,i(x)} 2 V(x) . 

xSB„ x,y£Z d 

Lemma 3.3. Fix a vector 9 in (R d ) n , £ > 2, denote by C\ the constant introduced 
in Lemma Vd .'A and let ip — ipg. There exists a finite constant ag, depending only on 
the transition probability p(-) , such that 

/log/^d/i n < - J fTp\ogiPdfi n - \og\A e \ n + 4(7i^(V?) + R{Q)Z 2 
for every density f : B n — > K with respect to ip d\x n . 

The estimates (|3.6|l permit also to prove the following bound. Recall that ft is 
the solution of the forward equation l|3.2[) and that /x„ is the counting measure on 

Bn- 

Lemma 3.4. Fix a smooth increasing function p : K + — > (1, oo) and a smooth func- 
tion A = (Ai, . .., An) : R+ (R d ) n . Let ^ t (x) = exp{A(i) • x} and let h t = ft/tpt, 
u t = hf^ 1 ^ 2 . There exists a finite constant ao, depending only on the transition 
probability p{-), such that 

J hf t} ip t dfin < J u\ log u 2 ipt dn n - (p(t) - 1) J u\ ip t d/j, n 

(P(*)-1), 



Pit) 



-S^(ttt) + R(X(t))p(t) / ui ip t dfin 



The proof of Lemma l3~4l relies on the estimates (|3.6|1 and follows closely the proof 
of Lemma 5.1 in [2]. 

We are now in a position to prove Theorem 13.11 Recall that ft is the solution 
of the forward equation Q3.2fl . Fix T > large, set q = 1 + (logT) _1 , g' = logT 
and consider a smooth increasing function p: [0, T] — > [9,9'] such that p(0) = q, 
p(T) — q' . At the end of the proof, pit) will be taken as a rescaling of the function 
[1 - (s/T)"]- 1 for some < a < 1/2. 

Following Davies g|, fix 9 = (9 1 , . . . , 9 n ) in (R d ) n , define ip t - B n -> R + by 

^(x)=ex P {-|i I ..x} ) 
denote 6>i p(i)/[p(*) - 1] by A»(i) and let h t = ft/ipt- 
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For a function g: B n — > R and 1 < p < oo, denote by 1 1 1 1 v =y> the L p norm of g 
with respect to the measure tp d/j, n '■ 

h%, P = E iffWi^w • 

A straightforward computation gives that 
|log|MU,.(t) = -^og\\h t \U t , pit) + ^_^-|||^|| f) p(t) . (3.8) 

Denote by and Ut/HutlH by v 2 . By Lemma [3.4I the second term on the 
right hand side of last formula is bounded above by 

A J v 2 \ogv^ t d^ n + ^\ og \\h t \\^ p{t) + R(t) (3.9) 
P(t) - 1 ^ f \ P(t) ~ 1 f 2 ; , 

77T2-2D^ t («t) 7^— / «i Vt^n , 

P(t) P\t) J 

where R(t) = R(\(t)). Notice that the second term in this expression cancels with 
the first term in the previous formula and that v 2 is a density with respect to the 
measure ipt d/j, n . By Lemma 13.31 the first term of this formula is bounded by 

m 



-J v^ t logVt^„ - log|A,r + AdfV^ivt) + i 2 R{t)) (3.10) 



for all I > 2. 

By definition of ip t , 



P(*)-l ; Pit) , , , 



P (t) ^ p(t) 2 

so that the first term of formula (|3.1U|) cancels with the fifth term of formula l|3.9|l . 
Denote by [a] the integer part of a real a. If we set I — l(t) as 



£(t) = 



I Pit) ~ 1 



a straightforward computation shows that the Dirichlet form in formula (|3.9|l can- 
cels with the Dirichlet form appearing in (|3.10(l . The inequality £(£) > 2 imposes 
conditions on p(t) that will need to be checked when defining p(t). 
Up to this point we proved that 

|log|NU, P ( t ) < -^log|A,r + 2i2(*) 
because £(t) 2 p(t) < p{t) 2 by definition of i(t). Integrating in time, we obtain that 
\\h T \U T , PT < ||/i |Uo,Poexp{-(nd/2)^ dt Ai og |fci + 2 J dtR(t)} . 

because £{t) 2 > [p(t) — l]/%C\p(t). By definition of the density /, ||/iolU ,po = 
f(z)ipo(z) 1 ~ p °/ po — exp{6* ■ z}. On the other hand, ||7it||v>t,pt i s bounded below by 
/t(x)V't(x) 1 ~ Pt ^ Pt = exp{# • x} for every x in B n - Moreover, since —p(t)/p(t) 2 = 
(1/P(*))', 



■logiSd) £dt A < log(8d 



SYMMETRIC EXCLUSION PROCESS 



15 



because p(0) — 1 + (logT) 1 , p(l) = logT. Finally, since p(t) is an increasing 
function, p(t)/[p(t) — 1] < 1+logT < 2 logT for T > e. Therefore, if we assume that 
< B/2\ogT for some finite constant B, R(t) < C(a , B)\\9\\ 2 p(t) 2 / (p(t) - l) 2 , 
where C(ao,B) — 2agM(ao-B) and 

cosh w — 1 
M(r) := sup 5 < coshr . 

Putting together all previous estimates, we obtain that 

/t(x) < C(n, d) exp{-0 ■ (x - z)} exp { - (nd/2) J dtA \ og M-l^ x 

provided |^-| < B/2\ogT. 

It remains to choose an appropriate increasing smooth function p : [0, T] — > [g, g'] 
which connects 1 + (logT) -1 to logT to conclude the proof of the theorem. Let 
q(s) = p{sT)/p{sT) - 1 and notice that q(0) = logT + 1, q(l) = logT/ logT - 1. 
With this notation, a change of variables and an elementary computation shows 
that the two previous integrals become 

(nd/2) fds^lo^-^zl _ (^/2)logT(J^--i-) 
Jo Q( s ) -Q{ s ) VlogT+1 logT/ 

and C(a ,B)\\e\\ 2 T dsq(s) 2 . 

Jo 

The second term of the first line is bounded by logT~ nd / 2 + C(n,d), which is 
responsible for the diagonal estimate of the density. 

Let g(s) — s~ a for some < a < 1/2. It is easy to show 

- 1 -ff'(s), g(s)-l f 1 , , , 2 

Defining g(s) = g(a + (6 — a)s) for appropriate constants a, b, we deduce that 
M x ) ^ ^^exp{-^.( X - Z )+C7(a 0)J B)||e|| 2 T} 

provided | < B/2 logT. An elementary computation shows that with this choice 
l(t) > 2 for all < t < T provided T is chosen large enough: T > C 2 (n, d). 

Fixx, let y = x-z and choose 6 = B(2 logT) _1 y/l|y| so that \9 hj \ < B/2 logT. 
With this choice, the expression inside braces in the previous formula becomes 
bounded by 

||x — z|| Ti C(a ,B)B 2 T 
2 logT + (2 logT) 2 
for every B > 0. Recall the definition of C(ao,B), change variables as B' = ciqB 
and minimize over B' to obtain that the previous expression is bounded above by 

a T ^||x-z||logT 



2(logT) 2 V a 2 T 
where $ is the convex conjugate of w 2 cosh w. This concludes the proof of Theorem 
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Proof of Theorem ll.il Th.eorcm. il . ll follows from Theorem ^, ll since the evolution 
of n random walks evolving with exclusion can be obtained from the evolution of 
n labeled random walks by just ignoring the labels. 
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